X 



DESY 06-043 |hep-ph/0604T60 

April 2006 



QCD Corrections to Semi-Inclusive Hadron 
Production in Electron-Positron Annihilation 

at Two Loops 



A. Mitov and S. Moch 



O 

o 

u _ 

< 



> 

O 

Deutsches Elektronensynchrotron DESY 
O i Platanenallee 6, D-15738 Zeuthen, Germany 

O 

iu: 



Abstract 

We report on the analytic calculation of the second-order QCD corrections to semi-inclusive hadron 
production in electron-positron annihilation. The calculation of the single-particle inclusive cross- 
section in time-like kinematics is performed in Mellin N-space and uses an algorithmic evaluation 
of inclusive phase-space integrals, based on the unitarity cutting rules and integration-by-parts. We 
obtain splitting functions and coefficient functions up to second order in the strong coupling ttg. 
Our results are in agreement with earlier calculations in the literature by Rijken and van Neerven. 



1 Introduction 

The direct observation of hadron production in e+e^ -annihilation offers unique ways to test predic- 
tions of QCD and has been measured in the past especially by the LEP experiments, see Ref. [1,2]. 
These data allow in particular for studies of single hadron production, fragmentation functions and 
their scaling violations. They also offer the possibility for determinations of the strong coupling 
constant tts at the scale Mz- 

In parallel to the deep-inelastic scattering (DIS) of leptons on hadrons, the fragmentation of 
the partons produced in the hard scattering process e+e^ ^ q,q,g into hadrons depends on the 
hadron's scaled momentum x. Perturbative QCD predicts the scale dependence and evolution 
of the scattering process, making higher order corrections mandatory for precision analyses. In 
the past, the next-to-leading order (NLO) results have been obtained in Refs. [3-6], while the 
coefficient functions at two loops, necessary for a next-to-next-to-leading order (NNLO) analysis 
were calculated [7-9]. In order to complete this program, the three-loop QCD predictions for 
the splitting function in time-like kinematics are required. The splitting functions governing the 
NNLO evolution of flavor non-singlet fragmentation distributions at third-order have recently been 
reported [10], while the corresponding singlet quantities are still missing. 

Our motivation for the present paper is at least two-fold. First of all, given that e+e -annihilation 
is of particular interest for the physics analyses of LEP data or a future ILC, we would like to pro- 
vide an independent cross-check on the computation of Rijken and van Neerven [7-9]. Secondly, 
we perform a first application of a new innovative method [11] to calculate higher order QCD cor- 
rections to single-particle inclusive observables directly in Mellin A'^-space. It is worth emphasiz- 
ing that the present method extends well beyond the framework of the operator product expansion 
(OPE) used in DIS [12-16]. This latter point is, in a more general perspective, rather important as 
it will allow for the efficient calculation of QCD corrections to many single scale observables. This 
can be done for Mellin moments either at fixed A'^ analogous to DIS calculations of Refs. [17-20] 
or by means of well established summation techniques [21] when full control over the analytic N 
dependence is kept. 

Finally, the present paper provides the means to thoroughly study the relation between the 
time-like QCD corrections to inclusive hadron production in e+e annihilation and their space- 
like counterparts, the DIS structure functions. At the leading order (LO) the Gribov-Lipatov rela- 
tion [22] suggests simple relations between the splitting functions in both kinematics, which do not 
hold beyond LO [3,23], see also Refs. [24-26]. Yet, the space- and time-like cases are related by 
an analytic continuation in x which has recently been used to obtain the NNLO flavor non-singlet 
time-like splitting functions [10]. In the present paper, we have also calculated in J = 4 — 2£ di- 
mensions so-far unknown terms at higher orders in e for the corresponding coefficient functions 
in e+e^-annihilation. These terms were subsequently used to check the analytical continuation 
between processes with space- and time-like kinematics [10]. 

The outline of the article is as follows. In Section |2l we set up the notation and the formalism 
for calculating QCD corrections to fragmentation functions. We also discuss the definition of the 
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(time-like) splitting functions and respective coefficient functions through order a^ . Section El 
briefly explains the method to calculate Feynman diagrams in Mellin A''-space for processes with 
time-like kinematics, with all details on the necessary master integrals and reductions given in 
Appendix A. In Section |4| we present our results for the fragmentation functions Fj^Fi and Fa 



through order a^ and with details on the mass factorization to extract the splitting and coefficient 
functions. The lengthy full expressions are deferred to Appendix B (//-space) and Appendix C 
(jc-space). We summarize in Section|5l 



2 The setup 

The subject of our calculation is the single-particle inclusive e+e^-annihilation, i.e. the process 

e+{h) + e-{k2) -^ V{q) -> H{ph) + X , (1) 

where V is a vector boson, i.e. V = y,Z, and {ki +k2f' = q^ = ^ > h its (time-like) four- 
momentum squared. The observed hadron with momentum pn is denoted by H and X stands for 
any hadronic final states allowed by quantum number conservation. 

Our goal is the derivation of the normalized double differential cross-section for the reaction in 
Eq. dU), 

— Z,"^ ^ =^(l+cos29)Fr(x)-F^sin29Fz.(x)-F^cos9FA(x), (2) 

Otot"^«cos9 8 4 4 

in perturbation theory, and including the quantum corrections to the fragmentation functions Fj^Fi 
and Fa through order (ttg ) in the strong coupling. In Eq. © o^ denotes the cross-section for 
producing the hadron H, x is its scaled momentum fraction: 

X = ^2 ' < ;c < 1 , (3) 

Q 

and 9 denotes the polar angle between the hadron and electron beam directions. The (total) frag- 
mentation functions Ft.Fi and Fa in Eq. Q originate from the transverse or longitudinal polar- 
ization states of the intermediate vector boson (Y,Z) or from parity violation of the electroweak 
interaction. Upon integration over 9 and x the asymmetric contribution Fa cancels and one arrives 
at the total cross-section integral 

-/ dxx-— —— = -— + -— = 1, (4) 

2./o Otot dx Otot Otot 

which represents the energy-momentum sum-rule. 

The evaluation of the fragmentation functions Fj.Fi and Fa in perturbative QCD is based on 
factorization. Up to power corrections suppressed by the hard scale Q, one can write the de- 
sired hadron level observables F as a convolution of collinearly renormalized, parton level cross- 



sections 'J with non-perturbative fragmentation distributions D. The explicit form of these rela- 
tions is: 



FiiA = oZliQ^) fi,q®D\_tj + Ji^g®Dg^H +^ai^(e2)-r;,„s®£>: 



«/ 



r(0)//i2w ^niis.- 



'h^Hi 



h=\ 



h=i 

where I = T,L and the symbol ® stands for the convolution integral, 

[f®g]iz) = fdx f'dyf{x)g{y)h{z-xy). (6) 

Jo JO 

The factors OjJ and A)^ are respectively the LO total cross-section and asymmetry factor for 
e'^e^ ^ hadrons. Moreover, 

of;M) = t^f\Q'). (7) 

h=\ 

The explicit expressions for aj^ and Aj^ can be found in Ref. [9]. 

The distributions D for a non-singlet quark (of flavor h), a singlet quark {q) or gluon g frag- 
menting to the observed hadron H are non-perturbative objects that are extracted from experimen- 
tal data at low scales Q of order 1 GeV and evolved to high scales by means of the time-like 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution equations [22,27-29]. 

As discussed above, the time-like splitting functions are presently known to NLO accuracy 
in the flavor singlet and to NNLO in the non-singlet case [10]. Thus far, this has permitted the 
determination of the fragmentation distributions with next-to-leading logarithmic (NLL) accuracy. 
Examples of sets of LL and NLL accuracy are [30-33]. In case the hadron H contains a heavy 
flavor, e.g. Z7-quark, the fragmentation distributions develop an additional perturbative component, 
a so-called 'perturbative fragmentation function' for the heavy quark [34]. That function allows 
the resummation of large logarithms of the ratio of the quark mass and the hard scale; it is currently 
known to NNLO [35, 36] which allows for the extraction of heavy quark fragmentation distribu- 
tions with NNLL accuracy provided all three-loop time-like splitting functions are available. 

The relation between the fragmentation distributions introduced above and the respective dis- 
tributions for specific flavors are. 



ns,+ 


= D,^^+Dj^^^-D%, 


ns,- _ 


= Dh^H-Dj^^fj, 


s 





(8) 

where q denotes any generic 'quark' flavor, while h stands for a specific quark flavor. The non- 
singlet distributions evolve with the corresponding combinations of splitting functions [10] (see 
also Eqs. (30)-(34) in Ref. [13]). 



Next we turn our attention to the evaluation of the hard partonic cross-sections J appearing in 
Eq. (|5l); the explicit results for these functions are discussed in Appendices B and C. The coefficient 
functions f are related to the partonic equivalents J of the corresponding observables F . For the 
construction of the 'bare' functions J one replaces the observed hadron H with an on-shell mass- 
less parton that can be q, q or g. Clearly, due to the incomplete inclusiveness of the observable dU), 
the partonic cross-sections J will contain additional coUinear singularities, i.e. these functions are 
divergent even after the usual UV renormalization. To that end one performs a so-called 'coUinear 
renormalization' or 'mass factorization'. One factorizes the coUinear singularities and effectively 
absorbs them into the fragmentation distributions D. As an intermediate regularization we work in 
d = A — 2£. dimensions [37-40]. In z-space and in the MS scheme [41], these relation take the form 
(with/ = r,L and J = T,L,A): 



0,ns(e,z) = 


= ify,ns(z)®r^^(e,z), 


^/,s(e,z) = 


= j/,s(z) ®r^^(£,z) +n/j/,^(z) ®r^^(e,z: 


^/,^(e,z) = 


= 27,,q{z)®Tgq{z,z) + Ji,g{z)®Tgg{z,z) 



(9) 

and in Section HI we explicitly present the corresponding expression in //-space. The form of the 
coUinear counter-terms T is universal, i.e. they contain only the time-like splitting functions. In 
the MS scheme, and in terms of the bare strong coupling, their explicit expressions can be found 
in Eqs. (4.25)-(4.30) of Ref. [9] and we do not repeat them here^ 

The partonic scaling variable z appearing in Eq. Q is the normalized energy fraction of the 
observed parton defined with respect to the total four-momentum q of the e+e system: 

z=^, 0<z<l. (10) 

Q 

This definition of z is consistent with the requirement that all partons are massless (see below). 

J in Eq. © above are the (finite) dimensionless partonic cross-sections with the coUinear (or 
mass) singularities factorized in the MS scheme. We have suppressed the explicit dependence 
on the factorization scale ^f\ throughout this paper we set it equal to the hard scale /i^ = Q^. 
Moreover, all functions appearing in Eq. Q are expressed in terms of the renormalized strong 
coupling. The relation between the bare ttg and the renormalized <x^{^r) couplings reads: 

The factor S^ = exp(£{ln(47r) —Ye}), where Ye denotes the Euler-Mascheroni constant, is an ar- 
tifact of dimensional regularization [37-40] kept out of the coefficient functions and anomalous 
dimensions in the MS scheme. Also, we have 
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'We would only like to caution the reader about the different notations for the definition of 8 and the normalization 
of the splitting functions used in these papers. In addition, the conventions for these counter-terms are the transposed 
of what is used in the more current literature. 



which is the 0{al) coefficient of the QCD P-function, Ca = 3, Tr = 1/2 are the QCD color 
factors, nf denotes the number of light fermion flavors and /jr stands for the renormalization scale. 
To simplify our expressions we will work with ijj^ = Q^. If desired, the explicit dependence on the 
arbitrary renormalization scale can be easily restored. 

To completely specify our observables, we have to clarify how the bare partonic cross-sections 
are defined (and calculated). This is done next in Section |3l 



3 The method 

3.1 General considerations 

Similarly to the DIS case, the bare partonic differential distribution J for the process Eq. ([T]) can 
be written as a product of a leptonic and hadronic tensors. The hadronic tensor Wt,v is proportional 
to the amplitude squared of the decaying vector boson, and depends in particular on the coupling of 
the latter to the fermion current. Throughout this paper we will only consider the case of massless 
quarks". In this case, as was detailed in Ref. [9], the calculation of the inclusive QCD corrections 
to the corresponding coefficient functions is independent of the vector boson being vector (V) or 
axial-vector (A) type. For its evaluation it is therefore sufficient to consider the decay of a vector 
boson V that couples to the fermion current as y"(1 — y^). 

The bare hadronic tensor W^v depends on two momenta: the one of the decaying vector boson 
q and the momentum of the observed parton p. Following the usual arguments for Lorentz and 
CP invariance one can show that W^v can be parameterized precisely with the three fragmentation 
functions 9^t,l,a introduced in Eq. ®. For the V-A couplings mentioned above (see [9] for the 
general case), these three functions can be expressed as projections of the hadronic tensor W^v 
(with J = 4-2e): 



a —2 \ q^ '^ p-q 



^r(z,e) = -^(-2^W;,^-^p>Xvi, 



p-q 

Our goal in this paper is to calculate the fragmentation functions Jt,l,a including the coefficient 
functions of order a^ . In fact, we even obtain the terms up to a^E^ in the mass-factorization, some 
of which we have used to check the analytical continuation of the space-like DIS results to the 
time-like region [10]. Terms of order a^E^ would also be needed in a future evaluation of the order 
a^ corrections to 9^t.l,a (see also Eqs. (I26b - (l39t for more details on that point). 

^Mass effects are known analytically through NLO [6] and, based on a numerical study, through NNLO [42,43]. 



Next, we explain how to construct the contributions of order a^ to the hadronic tensor W/jy. 
This tensor contains the contributions from all diagrams for the processes, 

V{q)^h{p)+hi{pi) + ... + h,{pc), c=l...k+l, (14) 

of the decay of the vector boson V with momentum ^ to a set of particles h,h\,. . .,hc. The dif- 
ferent final states (labeled by the number c of unobserved particles in that state) represent the 
contributions from the various physical cuts to the single-particle inclusive observable. After the 
Feynman diagrams contributing to each particular physical cut have been constructed and appro- 
priately added, one has to perform the required virtual and/or (real) phase space integrations. 

The approach we pursue for the evaluation of the required Feynman integrals differs signifi- 
cantly from the calculation in [7-9] . Our approach is based on the application of algebraic rela- 
tions based on integration by parts (IBP) [37,44-46] to cross-sections and it is performed directly 
in Mellin A^- space [11]. 

Our starting point is Eq. (IT3t : in the following we will use J to denote any one of the three 
bare cross-sections ^t,l,a- Each one of the functions iF has a perturbative expansion in terms of 
the bare^ strong coupling, 

^(^^^) = L&) f^'\n.t). (15) 

Throughout this paper, the letter n will be reserved for the Mellin variable of any function of z 
defined through 

/(„)= f dzff{z). (16) 

with n > 0. In particular, the total integral of a function corresponds to « = 0. This definition of the 
Mellin variable is the most natural choice for the calculational procedures detailed in the following. 
We will, however, present the final results for the corresponding finite partonic cross-sections J in 
terms of the conventional Mellin variable A'^ defined through 

ir(A^,e)= / Jzz^-V(z,e). (17) 

./o 

In view of the additional factor of z in the partonic equivalent of Eq. @ the relation between the 
two variables is 

n = N-2. (18) 

Each function J^^'{n,€) with k>Q contains contributions from a number of terms, corre- 
sponding to the different physical cuts of the process V — > /? -\-X at order a^: 

:FW(n,£) = £j([5^(n,£). (19) 

c=l 
^We find it more convenient to present results in terms of the bare coupling instead of the renormalized one, the 
reason being that we work with massless on-shell partons. Thus we only need UV renormalization, i.e. the one related 
to the coupling Eq. dlH . 



The functions y^}J{n,e) contain the full contributions from the process V ^>- h + X at order 
a^ where the inclusive final state X contains c unresolved partons. As described in Ref. [11], the 
construction of the functions !F}J{n,e) consists of the following steps: 

1. One constructs all contributing Feynman diagrams. The integrations over the virtual mo- 
menta FT ' are assumed implicit in the diagrams. 

2. One constructs the corresponding contribution to the tensor W^y by adding all relevant am- 
plitudes squared, and with the appropriate symmetry factors included. 

3. The above result is contracted with the appropriate tensor constructed from the J-dimensional 
metric and the momenta p and q as follows from Eq. (lOb . 

4. The Lorentz scalar constructed this way is integrated over the full phase-space of the c + 1 
partons (i.e. one also integrates over the full phase-space of the 'observed' parton). The 
measure for this integration is: 

d^ = (2p.qri2nYd{q-p-pi-...-p.)fl-^^d{pj). (20) 

We have defined pQ = p, and we have set Q^ = I for simplicity. The exact dependence on Q^ 
can be easily restored on dimensional grounds. The origin of the 'Mellin propagator' 2p ■ q 
is explained at the end of this Subsection. 

5. One applies the IBP identities to reduce each term (generally containing integrations over 
both real and virtual momenta) to a linear combination of a small number of independent 
'master' integrals. As a rule, at order a^ and for each particular cut c, one needs to construct 
and solve more than one IBP reduction; the number of the required reductions corresponds 
to the number of independent topologies for each set {k, c) . The 5- functions from the real- 
phase space are dealt with along the lines of Ref. [47], while the Mellin propagator is treated 
along the lines of Ref. [11]. 

6. Each master integral is a function of the Mellin variable n. Its n-dependence can be com- 
pletely extracted with the help of the difference equation the masters satisfy. The difference 
equations are obtained from the solutions to the IBP reduction (see also [16,48] for related 
discussions in the DIS case). 

7. One has to supply appropriate initial conditions for specifying the solutions of the difference 
equations. The most suitable choice is to evaluate the value of the masters at n = 0. This 
choice corresponds to the total integral of each master over z- Therefore the initial conditions 
are pure, £-dependent numbers. 

8. Following [11], we 'partial fraction' by performing an additional summation over n of the 
terms containing a propagator of the type ~ l/{l — 2p-q). This propagator is not linearly 



independent from the 'Mellin propagator' 2p-q as. it merely shifts the effective n in complete 
analogy to the DIS case, see e.g. Refs. [13,48]. 

For the evaluation of the transverse and the longitudinal functions it is sufficient to take the matrix 
Y^ as anti-commuting in J-dimensions. Special care is, however, needed for the evaluation of the 
asymmetric contribution JA.ns- 

For the evaluation of !FA,ns we follow the prescription of Larin [49]. The details about the 
implementation can be found e.g. in [13, 14]. In short, there are two important features: First, one 
replaces the axial-vector coupling y"y^ with the J-dimensional completely antisymmetric tensor 
g/^/poTy^y^y^ ^^^ l-]^gj^ ^gg^ j|-g coHtractlon properties with the second £-tensor appearing in Eq. (IT3t 
to reduce it to combinations of the J-dimensional metric tensor. Second, one has to multiply 
the resulting expression with additional renormalization constants which restore the axial Ward 
identity in the MS-scheme. These constants have an expansion in the renormalized coupling tts 
and in powers of £. They have been computed to three-loops in [50] and take the following form: 



Z, ^ 1 + ^^c^ I _4-10£-F (-22 4-2^2)82-^0 (£3)} 






+ (^(132-48C3)C/+(-^+48C3)QC^ + ^n^Q^£+0(£2)} + o(a3) 

Za = l + (^)'l{fQC.J.,Q} + 0(a3). (21) 

We would like to conclude this Subsection with a comment on the origin of the 'Mellin propa- 
gator' Pm = 2p-q appearing in Eq. (OUb . As was detailed in [1 1], to construct the bare distribution 
[F (z, £) in z-space, one has to integrate over the full phase space of all final states particle and in- 
sert the additional factor 5(z — 2p-q). If one Mellin-transforms this expression before the required 
phase-space and virtual integrations are performed, one gets schematically (see also Eqs. (ITOb . 

9{n,e) = (...) X / dzz"d{z-2p-q) 
Jo 

= (...) X {2p-q)\ (22) 

i.e. the factor of the 'Mellin propagator' raised to a symbolic power n that appears in Eq. (l20b . The 
factor (. . .) stands for the various propagators (including possibly additional powers of Pm), the 
measures for the real and/or virtual integrations, etc., but contains no dependence on z or n. This 
procedure is completely analogous but more general than the corresponding DIS case [12-16], 
which relies on the OPE and the method of projection to directly expand propagators in powers of 
{2p ■ q/q^)". There Eq. (|22ll is effectively realized by mapping any Feynman diagram to Mellin 
moments with the help of a suitable projection operator [51,52]. 

In the following we will present the specifics of the implementation of the above procedure for 
the evaluation of the contributions at orders a^ and a^ . 
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3.2 Order o (as) 

Here we discuss both the derivation and the results for all independent contributions at order (tts) . 
This order is the NLO result for the transverse and the asymmetric contributions, but the LO result 
for the longitudinal function, since it vanishes at order (oc^). 

At order 0(as) in the perturbative expansion, one has to consider only two cuts: one corre- 
sponding to real gluon emission (where the final state is (qqg)) and the one with the virtual correc- 
tion to the Vqq vertex (with the final state being (qq)). The evaluation of these two contributions 
is fairly different. We discuss first the contribution from the virtual corrections. 

From the kinematics of the tree-level decay process V ^ qq it is clear that the contribution of 
the purely virtual corrections to the functions J are of the type const {e)d{l — z) in z-space which 
corresponds to a n-independent constant in «-space. Therefore, the contribution from this cut is 
completely determined by (twice) the real part of the time-like form-factor [53-55]. 

All non-trivial n-dependence at that order comes from the real-emission diagrams. The corre- 
sponding diagrams can be found in Ref. [9] and we do not repeat them here. As was outlined in the 
previous Subsection, for the evaluation of all contributions {T,L,A) for both a quark and a gluon, 
one needs to construct a single topology consisting of five 'propagators' (see also [11]). One of 
the arguments, of course, corresponds to the Mellin propagator 2p ■ q. We performed the required 
IBP reductions and obtained a single n-dependent master integral, which can be found e.g. in [1 1]. 

To perform the reductions resulting from the IBP identities, we have used the program AIR [56] 
which is an implementation of the so-called Laporta algorithm [57-60] in MAPLE. AIR contains 
also a routine which conveniently and automatically maps the constructed diagrams into the master 
integrals of the reductions (a single master at this order). Following this simple procedure, one can 
map the whole problem of the evaluation of any one of the functions J^ at order a^ to the single 
master multiplied by a rational function of n and £. As discussed above there is another type of 
contributions containing a propagator that is not linearly independent of the Mellin propagator 
Pm = 2p- q, i.e. contributions of the form: 



pn oo 

(■■■)t^=L(---)Pm-"' (23) 

k=0 



1-Pm 



where dots stand for any other, linearly independent propagator and the integration measure over 
the real momenta. Since the term on the right hand side of Eq. (E^ is of the usual form, the results 
from the solutions to the IBP reduction can be applied (with n replaced hy n-\-k) and then summed 
over k. To illustrate that point, we present the corresponding term of the order a^ contribution to 



the function J., 



T. 



Cf £ z{l-z){A£^ -A£^k-5A£^ + 3\zk + lAz + zk^ -Ak^ -2Qk-2A) 
X n'-er-r(t+l-2e) (2^, 



r(2-2£)r(4-3e + yt) 



which can be easily summed up in terms of F- functions. This way, one arrives at a very compact 
result for the bare {T,L,A) partonic cross-sections at order a^ valid to all orders in £. The explicit 
results, expanded to sufficient powers in £, can be found in Section|4| 



3.3 Order 0{a: 



At this order one has to consider three different cuts: double- virtual corrections with final state 
(q^q), one-loop virtual corrections to the final state (q,q,g) and, finally, the cut with double real 
emission. The latter consists of the following final states: {q,q,g,g); {q,q,q,q) and {q,q,q',q'). 
Depending on the gauge choice for the polarization of the external gluons one may also have 
to consider external ghosts. Again the corresponding diagrams can be found in Ref. [9] and we 
are only interested in the evaluation of the diagrams with real emissions since the diagrams with 
two-loop virtual corrections produce only constant terms in Mellin n-space. These purely virtual 
contributions can be obtained from the one- and two-loop time-like form-factor [53-55]. 

To cover all possible diagrams for the evaluation of both quark and gluon production we con- 
struct seven topologies for the double real emission cut, and five topologies for the real-virtual cut. 
After symmetry considerations, we arrive at a total of six real-real and five different real-virtual 
masters. These n-dependent masters satisfy difference equations that can be read off the com- 
pleted reductions. For completeness, we have presented both the definitions of the masters and the 
difference equations they satisfy in Appendix A. 

As can be seen there, the structure of the equations is simple. The simplest masters decouple 
and satisfy homogeneous equations, while the more complicated ones satisfy first order difference 
equations with the non-homogeneous part comprised by simpler, explicitly known masters. Such 
first order equations can be easily solved in closed form to all orders in n and e, see e.g. Ref. [13]. 
If we pursue this approach, the most complicated solutions we encounter are Appel functions of 
unit arguments or hypergeometric functions of unit argument. 

However, we decided to follow a different path. One can also obtain the solutions of the 
difference equations after an expansion in powers of £ using the methods of symbolic summation 
and the packages Summer [61] and XSummer [62] in Form [63]. Then it is very easy to 
solve the equations this way given that previously even three-loop master integrals have been 
computed [14, 55, 62]. In this approach, one is required to supply the initial conditions (for n = 0) 
beforehand. This is to be contrasted with the all-order in £ calculation where the initial condition 
factorizes completely. 

Only one master, R2{n) deserves special consideration. It formally satisfies a second order 
difference equation as can be seen from Eq. (IA.13I) . A closer inspection, however, reveals that it 
is a second order difference equation of defined parity (thus a n— I term is absent). Therefore 
one can write this equation as a first order difference equation for a 'new variable' n' = n/2. In 
doing so we could solve the resulting difference equation in terms of vFg-type functions of unit 
argument. These functions contain half-integer (and n/2) parameters and are not simple to expand 
in £. If one, however, solves this master as an expansion in powers of £ no particular complications 
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arise besides the fact that one has to supply two initial conditions (for n = and for n = \) for any 
second order equation. 

Next we address the calculation of the initial conditions for the masters. Our first step is 
to perform new IBP reductions in each of the real-real and real- virtual topologies setting n = 
from the very beginning. From this fixed order reductions we obtain relations between the initial 
conditions for the masters integrals. In the end, we find that a total of seven are independent. Most 
of these integrals are very easy to compute directly to all orders in £. The most involved is the 
initial condition R(i{0) given in Eq. (IA.23I) . For its evaluation (up to weight 6 in values of the 
Riemann zeta-function) we have used the approach of Ref. [64] based on the optical theorem and 
the results of Ref. [65] for the higher order £ terms of the non-planar three-loop two-point function. 
With the help of these additional fixed-/? runs we can derive a relation (see Eq. (IA.23I) ) between 
the values of the real master i?2 at n = and n= 1 as mentioned above. 

The last point that deserves special attention are terms ~ 1/(1 — 2p ■ ^). 'Partial fractioning' 
of these leads to an additional sum as explained above, cf. Eq. (l23t . In performing the required 
symbolic summation we employ the following strategy. We explicitly separate sums from 1 to 
oo over ^-independent terms, as we systematically ignore constant terms from the evaluation of 
the diagrams with real emissions as well as the constant terms from the purely virtual corrections. 
After we complete our evaluation we can restore these constant terms from the requirement that the 
total cross-section is reproduced (see Eq. dUl). We wish to emphasize, though, that this procedure 
is simply done to economize on the necessary algebra and by no means represents any principal 
drawback of our approach. 

The above comment applies to the transverse and asymmetric partonic cross-sections. Since 
the longitudinal cross-sections do not contain any constant terms, this procedure uniquely fixes 
the missing constant contributions in the transverse functions. It also provides the coefficient of 
the 5-function contributions to the asymmetric functions. To that end, one uses the fact that the 
difference of the transverse and asymmetric functions when expressed in z-space, does not contain 
5- functions and singular -i-distributions. Equivalently, the difference of these functions vanishes 
in the 'soft' limit N ^ oo. 



4 Results 

Let us now present the results of the calculation. As explained above, mass factorization (or 
coUinear renormalization) predicts a specific structure for the final result, which we write out up 
to second order in the strong coupling a^ = as/(4;r). Following the conventions from Eq. (ITSl) . we 
present the general structure of the result expanded in e directly in Mellin A^-space. According to 
Eq. (UTJ), the Mellin moments of the splitting functions are defined by 



P(N) = Cdzz^-'P{z) = -Y(A^), 

JQ 



(25) 
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and we note the (conventional) sign for the relation of the splitting functions to the anomalous 
dimensions. In the following, all products are to be evaluated employing the algebra for harmonic 
sums [13,61]. 

The zeroth-order contributions, with jFy: being suitably normalized, read 

^(0) _ JO) _ J (0) _ (0) _ (0) _ (0) _ (0) _ 

Note that not all functions are independent; on general grounds one can show that [9]: 

^T,q = +^r,q, 
^A,q = -^A,q, 
!FA,g = 0. (27) 

For the applications of the present study we present the amplitudes at the first order in ttg up to 
terms of the order £^, yielding for iFr 

yO) ^ _?p(0)+40+e4'] + e^4'i, (29) 
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for the longitudinal jFl 



^S = cS+£4;j+£^4;j, (30) 

^(1) _ Jl) , pJl) , p2Jl) (21) 

-^L,g - ^L,g + ^^L,g + ^ '^L,g^ '<^'^) 



and for the asymmetric J^ 



^Z = 4<' + 4Ve4',', + e^4',',. (32) 

Correspondingly, the a^ contributions where the non-singlet and singlet quark amplitudes dif- 
fer for the first time, are required up to order £. These quantities are given by 

(2) _ 1 f viO) f niO) , aM 1 /d(1)+ loJI)d(O) 



+ c(2) _^(l)p(0) , J^(2) _^(l)p(0)l .33^ 

(2) _ J_rp(0)p(0)\_^rp(l)s (1) (0) 

-^r,ps — 2e2 1 'IS ^1 J 2£ I '''' "^'r^g^qg 

+ C^2) _l^(l)p(0) 1 12^(2) _^(l)p(0), .34. 

^ ^r,ps 2^'g IS ^ 2 I ^'P** '^r,gMg f ' v-J^; 



(2) _ lrp(0)/^p(0) (0) g\1 Ir (1) (1) (0)(1) ( 

-^r,g — e^ I sq \^ qq ^^gg ^^^Oyj ^ el ^^'-r.q'^gq ^'-r.g^g 



,(0) 
-t-Cj,g zaj-qZ-gq "7',g^gg ^^'^l^r.g ^'^r,q^gq '^T,g^gg f ' W^; 
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(2) _ _l/.(l)p(0)\ ^.(2) _.(l)p(0)^e/.(2) ^a(1)p(0)' 
L,ns p 

(2) _ 1 /J1)d(0)\ L J2) 1 (l)n(O) l,/o.(2) J1)d(0)' 



^(2) _ _£L(l)p(0)\ , ^(2) _£ (1) (0) 1 r (2) _ (1) (0)| 

^(2) _ 1 /oJl)p(O) , Jl)n(0)l , (2)_. (1) (0)_ (l)p(O) 
-^L,g — g ^^'■L.q^gq "T"'^L,g^gg J + '^L,g ^"L,q^gq "L,g^gg 



+ e^S-24>i?-4>£^r (38) 



-^A,ns — 2£2 I '''' V '•'' /J 2£ I +^M,q^' 



(1) p(0) 

A.q^qq 



+ .(2) _ (l)p(O) r (2) _Ai)piO)\ .30. 

^ ^A,ns "A,qMq ^ '^ \"A,ns ^A,qMq J • W^*; 

The results for the splitting and coefficient functions at order 0(as) and 0{a'^) are given both in 
A'^-space (Appendix B) and x-space (Appendix C). The precise definition of the various splitting 
functions can be found in Eqs. (30)-(34) in [13], see also [66] for a more detailed discussion on 
that point. 

Several comments are in order: First of all, as we have performed the calculation in Mellin 
space, all jc-dependence is recovered from the A'^-space results by an inverse Mellin transformation, 
which expresses these functions in terms of harmonic polylogarithms [67]. The inverse Mellin 
transformation exploits an isomorphism between the set of harmonic sums for even or odd A'^ and 
the set of harmonic polylogarithms (see also Appendix B). The algebraic procedure [13,67] is 
based on the fact that harmonic sums occur as coefficients of the Taylor expansion of harmonic 
polylogarithms. 

Our results for the finite terms in £ agree with the ones in [7-9]. We have found several 
misprints in these references and would like to take the opportunity to point out these typos in the 
original manuscript of Ref. [9] (employing the notation of the original reference). In Eq. (A.6) of 
Ref. [9] for Cj |^ there should be a replacement of the term 

C/{\+z){...-3\n^z...) -^ C/{l+z){...-3\nh...), (40) 

in Eq. (A.8) of Ref. [9] for ^^^^^^'"^ of the term 



(C/--QC^)(^ + ...llnz ^ (c/--QC^)(- + ...)lnz, (41) 



in Eq. (A. 10) of Ref. [9] for 4^^ of the term 
CpTf 



11 . 

+ —\nh+... — > CpTf 
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...+ l^Wz]+..., (42) 



and in Eq. (A. 15) of Ref. [9] for cJ^J(2),id ^^^^^ ^^^ 



L,q 

(C/-1qc^)(325i,2(1-z) + ...) — (C/-^QC^)(325i,2(-z) + ...)• (43) 



Finally, in Eq. (17) of Ref. [8] for C^ '"^ —Cj '"' one should replace the term 

C/(-^ + ..A\nz — > C/(-^ + ..A\nz. (44) 



Beyond the coefficient functions at order (a^ ) we have also obtained the terms of higher order 
in the £ expansion, specifically a'^^^Z?^^^ and a^^^. In the non-singlet case, we have found them to 
agree with the predictions based on the analytical continuation proposed in [10]. However, as these 
expressions (a^^', b^^' and a^^') are particularly lengthy and have no direct physical application, we 
refrain from writing them out explicitly here. 

For future use and for uniformity of the notations, we present in the Appendices B and C also 
the explicit expressions for the one- and two-loop time-like splitting functions. Our calculations 
agree with the known results [3, 4] (see also [66]). This statement, though, is subject to one 
qualification. As we are considering in Eq. ([T]) only the decay of a vector boson V, we have 
no access to the two-loop (time-like) splitting functions i^qg and Pgg with the set of Feynman 
diagrams considered. To do so, we would actually be required to compute also the decay of a 
(fictitious) classical scalar (|) that couples directly only to the gluon field via (|)G"yGfl^. In time-like 

kinematics, this approach has been used for instance to derive Pgg in [68]; see also Refs. [15,16,18] 
for the analogous considerations in the space-like case. An inclusion of the ^CLGa"^ coupling 

is straight forward and would allow for the determination of Pqg and Pgg (or rather its Mellin 
transform) to the desired two-loop accuracy. 



5 Summary 

We have calculated the O (a^) corrections to the transverse, longitudinal and asymmetric fragmen- 
tation functions for both quarks and gluons in semi-inclusive e+e^ -annihilation to hadrons. Our 
calculation confirms the results of Rijken and van Neerven [7-9] and we have taken the opportu- 
nity to correct several typographical errors in these papers (see also [10]). Our results constitute a 
strong check on [7-9], in particular since we have used a rather different technology and obtained 
them directly in Mellin A^-space following the proposal in [11]. Thus, our calculation represents 
the first example of a single-particle inclusive observable beyond the well established DIS frame- 
work [12-16] that is computed analytically in Mellin A'^-space. 

The coefficient functions presented in this paper contain the NNLO corrections to the transverse 
and asymmetric fragmentation functions Jj and Ja, and the NLO corrections to the longitudinal 
ones, Jl- After the complete singlet three-loop time-like splitting functions become available (the 
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non-singlet case has recently been reported [10]), one will be able to study light and heavy quark 
fragmentation at NNLO [34-36]. This important class of observables will allow precise extraction 
of fragmentation distributions from LEP data and is yet another motivation for the realization of 
the envisioned high-precision Giga-Z option of the future LLC. Due to the process independence 
of the fragmentation distributions, they can be further applied to other processes like hadro- and 
photo-production. 

In addition, the present paper provides the means to thoroughly study the relations between the 
time-like QCD corrections to inclusive hadron production in e+e annihilation and their space- 
like counterparts, the DIS structure functions. Our calculational approach easily allows us to obtain 
higher powers in £ of the bare partonic cross-sections at order (cx^ ) , which had not been computed 
before. Thus, we could provide important cross checks on the procedure of Ref. [10] based on an 
analytic continuation in x between observables with space- and time-like kinematics . Of course, the 
aforementioned higher terms in £ will also be needed for a future evaluation of the QCD corrections 
to e+e^-annihilation at order 0{al). Finally, one particularly appealing feature of the A^-space 
approach is the small number of master integrals that have to be evaluated. Moreover, with the 
boundary conditions in A^-space being kinematics independent, the corresponding integrals may 
also be of relevance in other circumstances. As a matter of fact, some of them had been considered 
before in a different context [64]. In the present paper we have extended these results to higher 
powers in £. 

Among the prospects for future developments and applications of our results and methods 
are explicit three-loop checks of the splitting and coefficient functions in semi-inclusive e+e^- 
annihilation, for instance by computations of fixed-A'^ Mellin moments. Also QCD corrections to 
many other single scale observables can be considered. 

Form files of our results can be obtained from the preprint server http : / /ar Xiv . orq| by 
downloading the source. Furthermore they are available from the authors upon request. 

Acknowledgments: We are grateful to A. Vogt for useful discussions. A.M. acknowledges 
support by the Alexander von Humboldt Foundation. The work of S.M. has been supported in part 
by the Helmholtz Gemeinschaft under contract VH-NG-105. 

Note added: The (un-)polarized coefficient functions up to two loops have recently been trans- 
formed to Mellin A'^- space in [70]. 

Appendix A: Master Integrals 

In this Appendix, we present the complete list of master integrals, the corresponding difference 
equations in Mellin n-space and the respective boundary conditions. We omit the discussion of 
the so-called purely virtual contributions, which are known since long. For the calculations of 
two-loop form factors in QCD we refer to Refs. [54, 55, 69]. 

Let us start with the master integrals Vi (n), . . . , V5(n) for the so-called real-virtual contribu- 
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tions. These masters can be defined through the following object: 

[i,...J]= (A.l) 

^(3Yee) /■ 1 

— ^ / d'kd'p,d'p2dipi)8ipl)d{{q-pi -P2f) {2q-{q-pi - Pi))" -^-^ . 

% J ri... ij 

where q-q = I and the propagators are Pi = k^, Pi = {q — Pi + k)^, Ps = {q — k)^i A = {pi + 
P2 + kf,Ps = {pi-kf,Pe = {p2 + kf,Pn = {q-p2f. 

The n-dependent real-virtual masters are defined as: Vi(n) = [1,2] , V2(n) = [1,3] , ¥3(77) = 
[1,4] , V4(n) = [1,2,5] , V5(n) = [1,2,5,6,7]. These masters satisfy the following difference equa- 
tions: 

(n + 2-4£)Vi(«)-(n + l-3£)Vi(R-l) =0, (A.2) 

(n + 2 - 3£)V2(n) - (n + 1 - 2£)V2(n - 1) = , (A.3) 

(n + 2-4£)V3(n)-(77 + l-2£)V3(n-l) =0, (A.4) 

(n+l-2£)V4(n)-nV4(n-l) = (A.5) 
-(l-3£)(l-2£)ivi(n-l) + (l-2£)2iv2(n-l), 

(n - 1 - 4£) V5 (n) - (n - 1 - 2£) V5 {n-\) = (A.6) 

(l-3£)(l-2£)(n + l-4£), . 9x 1 . ^ 

- (n-3£)(n-l-3£) ^ 27£^ - ll^B + 5£-n + .^)^VKn- 1) 

(l-2£)2(n-4£)(n + l-3£) 1 ^ , ^, , ^ ,^^ , ^, 
+ ^ n-2z ^^V2(n-l) + (n-2£)V4(/7-l). 

The boundary conditions Vi (0) , . . . , V5 (0) at n = read, 



Vi(0) = V3(0) = (A.7) 

V2(0) = (A.8) 

^^^^^^.(^^'^ ^^rUp2rl597 13 289 X 3/I2359 221 3111 
£8'^T6+H^-T6W+' l-6^-Y^3-^W+' l^2r"T^^^""6^^^ 

897 „ 2\ 4/88629 303„ 2379„ 27149„ 22r „ 15249^2 
+ 320^^ J +n^5r " ^^^ " ^r^3 " ^28"^^ + T^^^^^ + "640"^^ 

5/603871 515U 2076U 169^2 210103^ 3757^^ 16415U 2 

+ ' l^T^-^0"^^""6^^^ + T^^^ -^5^^^ + ^r^^^3 + ^280-^^ 
12949 3 

~ 4480 ^^ 
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V4(0) = (A.9) 

11 11 K f 11 5^ 49^ \ 2f 439 93^ 363^ 3^ ^ 

./ 2229 75^ 583^ 2141 13U ^ 1857^ 2\ a{ 10527 6303^ 

308U ,,^2 llllU 2317^^ 15147^2 421 3 

V5(0) = (A.10) 

15 15 1/5 133. \ 1/ 133. 193^ \ ^^ 133^ 5477^2 

193^ /^^ 4545^. 133^ 5477^2 193^ 2303^ \ 2/.^ 

--C3+e(20+^C2C3--C2 + ^C2^--C3-^C5)+BX40 

+ ^^2^3- 133^2 + ^0-^2 -^3440-^2 - 193^3 + ^^^3 -^^5 

which we have given in terms of the Riemann zeta- function consistently up to weight 6. Previously, 
the result for V5(0) has been obtained to weight 4 in Ref. [64]. 

Next we present the master integrals Ri(n),. ..,R6(n) from the so-called real-real contribu- 
tions. These masters can be defined through the following object: 



{/,...,;•} = —^Jd'pid'p2d'p3^{plWplMplM{q-pi-P2-P3?) (A.ll) 
X {2q -{q-pi-pi- P3)T — , 

where ^-^=1 and the propagators are 2i = (^-p3)2, Q2 = {q-P2f, Q3 = {q-Pi-Psf, 24 = 

{q-p2-P3f. Q5 = {q-P\-P2f. 26 = (Pl+P2)^ Ql = {P\+P3f. Qs = {P2+P3f- 

The n-dependent real-real masters are defined as: Ri(«) = {— }, R2(") = {I72}, R3(n) = 
{2,6}, R4{n) = {1,2,3,5}, R5(n) = {1,2,6,7}, Rein) = {3,4,7,8} and satisfy the following 
difference equations in n together with boundary conditions at n = 0: 

(n + 3-4£)Ri(n)-(n + l-2£)Ri(n-l) =0, (A.12) 

(n+l-2£)(n + 2-6£)R2H-(n-l)(n-4e)R2(n-2) = (A.13) 

^ ( 1 - £) ( 1 - 3e) (2 - 3e) (2n + 1 - 6£) , 



(n-3£)(n+l-3£) 



■Ri(n-2), 



(n + l-2£)R3W-.R3(.-l) = - ^^-^^^^^;f2^37^-^^^ lRKn-l), (A.14) 
(n-l-4£)R4(n) + (n-l-2£)R4(n-l) = (A.15) 

(n - 2e) (n - 3£) (n - 4£) (n - 1 - 3£) (n - 1 - 4£) (n + 1 - 3£) V 
- 1098£^ + 3282n2£4 - 279£4 + 1407n£4 + 52£3 + 269n£3 - 1 135n3£3 - 707nV 
-4lne^ + llAn^e^ + 5e^ + llln^z^ - 98nV + 16n3£ -2\n\ + 1 Inh -nz- 23n^£ 
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A .3 , ..5 , 6A 1 o ^„ ^, (n+l-2£)(n + 2-6£)(n-6£)l 



-n'-n' + n^+n'']^Ri{n-\) + ^ '-^ '-^ ^-R2(n) 

/ £2 ' ' n-At £ 

(n-2£)(n+l-6£)(n-l-6£)l , ^, 

; 'a K2(/7 — I) , 

n-\ -At £ ^ ^ 

(n-l-4£)R5(«)-(n-l-2£)R5(n-l) = (A.16) 

(l-3£)(2-3£)(n + 2-4£) 



('l20£^ - l5Am^ - 142£3 + llnh^ 



+ 104n£^ + 23£2 - 10n£ - 25n^£ - 14n^£ -z + r? + n^ + 2n^]^Kx{n-l) 



(n-3£)(n-4£)(n-l-4£)(n+l-3£) 

1 

£2 

^ NT. / ,N (n + l-2£)(n + 2-6£)(n-6£)l , , 

+ 2 n-2£ R3 n-1 +^ ^^ 1 ^^ —^2{n) 

n — Az £ 

(n-2£)(n+l-6£)(/7-l-6£) l 

n-l-4£ £ ^^ ^' 

(n-l-4£)R6(«)-(n-l-2£)R6(/7-l) =0. (A.17) 



Ri(0) = (A. 18) 

1 71 2/3115 7 ^ \ 3/109403 29^ 497^ \ 4/3386467 2059 



^V^^ I 20736 96^^ 384 ^V ^^ I 124416 576^ 



96 576 V3456 64 "V V 20736 96" 384 "V V 124416 576 
21805^ 291 ^ 2\ ./96885467 421 90335^ 765821 203^^ 



2304 " 1280'' / V 746496 160'' 3456 " 13824 " 64 
6887 2\ 6/2631913075 29891 3172687 841 2 _ 23705269 

^2560^^ )^^ \ 4478976 960"^^ ~ 20736 ^^^192^^ 82944 ^^ 

14413 60431 2 15089 3 
+ 384 ^^^^ ^ 3072 ^^ ~ 53760^^ 
R2(0) = (A. 19) 

IK / V 9^ 7. \ 2/ 119 63. 87. 97. 2\ 3/ 199 
-8 + 8^^ + <-4 + 8^^ + 8W+^ (-^ + ¥^^+16^^ + 80^0+^ (- — 

207^ 163^ 139. 21U . 679. 9\ ./ 4617 1449. 1585, 



45.2 3445. 1477.. 11931.2 1141. 3\ 



R3(0) = (A.20) 

11 11 ,/ 7^21 X 2/ 439 29 231 X 3/ 2229 319 

1617. 873. 2\ 4/ 10527 1263. 2233. 9219. 609.. 

+ Yr^^~320^' J+' \ ^ + ^0"^^ + ^~^^ + Yr^'"T^^^^^ 

9603^ 2\ sf 47389 13893^ 1273U 841 2 46809^ 6699^^ 



320 "/ V 8 ' 40 "-^ 8 ""^ 16 ""^ 16 " 16 
67221 2 15089 3 

320 ^^ ^ 4480 ^^ 
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=3 



R4(0) = R5(0) = (A.21) 

13 13 1/3 83 \ 1/3 185 83 \ 3 185 83 457 2 
£432^i3T6^£2U" 64^V + i U " ~yi^^ ~ 32^V + 2 " l6"^^ " 16^^ " 256^^ 
/ 15899^ 185^ 83^ 4137^^ 457^ 2\ if, 15899^ 185^ 

6667^ 2 83. 4137^ ^ 457^ 2 692333^ jX 

R6(0) = (A.22) 

15 15 1/5 105^ \ 1/ 141 105^ \ ^^ 141 105^ 

969^2 /o^ 1119^ 14K 105^ 2961.^ 969. 2\ 2/.^ 1119^ 

,.,. 4005.2 ,r.cr 2961^. 969^2 15121. 3\ 

Again we have given the results consistently up to weight 6 in the Riemann zeta-f unction. The 
difference equation for R2(n) which is of second order needs actually two boundary conditions. 
The value for R2(l) is obtained from an additional fixed-n reduction 

R2(l) = ^Ri(O) ■ (A.23) 

The values for R4(0), R5(0) and R6(0) have been obtained before up to weight 4 in Ref. [64]. 

Appendix B: The exact Mellin-space results 

Here we give the exact Mellin-A'' expressions for the time-like splitting functions and the coefficient 
functions cj, cl and ca up to second order in as = as/(47r), expressed in terms of harmonic sums, 
recursively defined by [61] 

S±,n,,m^....,,nM = £ ^ S,„,,...,„,, (0 , S{N) = 1 , (B.l) 

and we employ the notation 

N±5^ = Sm{N± 1) , N±iS^, = Sff,{N±i) . (B.2) 

The well known results for the LO anomalous dimensions (cf. Eq. (E51) ) in A/^-space are: 

yS^A^) =q|-3 + (2N_ + 2N+)Si| , (B.3) 

i^^{N) = «J2N_ + 8N+-4N+2-6l5i , (B.4) 



19 



yg) (AT) = q| _ 8N_ + 4N_2 - 2N+ + elsi , (B.5) 

Yg\A^)=Q|-y + (-8N_+4N_2-8N++4N+2 + 12)5i| + ^n^. (B.6) 

For the NLO anomalous dimensions (of. Eq. ^I5\i ) a distinction between even and odd Mellin 
moments is in order. We give the respective expressions for the even Mellin moments of Yns in 
Eq. dnnil and y^^^\ yjg^ y^q^ y^g^ in Eqs. (iB^- dBTTIb and for the odd Mellin moments of yll^" in 
Eq. (IB. 81) . Note also, that we do not obtain yqg^ and ygg since we only consider vector boson decays 
(see discussion in Section|H). Nevertheless, we also quote these quantities here for completeness. 

yi!J^{N) = (B.7) 

Cf^I ( 8- (16^2 + 4)N+ - (16^2 + 4)N_ )5i-(4N+ + 8)52- (4N_+28N+)53 + 



(-16N_ - 16N+)5i,_2 + (8N_ + 8N+)5i,2 + (8N_ + 8N+)52,i + 165_3 + 24^2 " ^ K 

r /302 112 302 \ / 22 22 \ 

C^q|-85_3+(-^N+- — + — N_j5i+(--N---N+j52 + 8N+53 

^ + (8N_ + 8N+)5i,_2} + 

44^, 44^^ 16\ 1 /4^^ 4, 



Cpnfi { -— N+-— N_ + — )5i + -+( -N_ + -N+ )52 



yi^^~{N) = (B.8) 

C/| 165-3 - ((16^2 + 36)N+ -72 + (16^2 + 36)N_)5i + (16N_ - 20N+ - 8)52 - 
(20N_ + 12N+)53 + (8N+ + 8N_)52,i - (16N_ + 16N+)5i,_2 + (8N+ + 8N_)5i,2 - 

r /446 208 446 \ / 46 2 \ 

C^q|-85_3+(-^N-- — + — N+J5i+( -yN- + -N+J52 + 8N-53- 

17 1 

— + (8N+ + 8N-)5i,-2[ + 

44 44 16\ 1 /4 4 

yN.-^N_ + - 5, + -+ -N, + -I 



Cpnfi 1 -— N+-— N_ + — )5i + -+( -N+ + -N_ )52 



li!^\N) = (B.9) 

208 _ 224^^ _ 80^, 80. 
_N+, + 48--N_2 + - 



Cp^fi ( -^^N--^^N+2 + 48-— N_2 + — N+)5i + 



-YN+-yN+2 + 16 + 20N_)52 + (8N_-8N+)53 



20 



i^JiN) 



(B.IO) 



C^Rj|(-86N++40N+2 + V0-24N_)5i + (10N_ + 8N+2-18)52+(l6N+2+12 



4N_-24N+ 53 + (16N+-20-8N+2+12N_)Si,i+ 96N+ + 24N_-48N+2- 



72 ] 5i,2 + (24 - 8N_ + 16N+2 - 32N+)52,i + (24 - 8N_ + I6N+2 - 32N+)5i,i j | + 

Qn^{ ((y-16C2)N_ -(64C2 + ^)n+ +(32^2-^) N+2 + 48^2- 

80 2224\ / 16 16 \ , 

-5-N_2 + -5- ) "^1 + ( - yN+2 + 40 - 40N+ + — N- j 52 + (8N_ - 48N+ + 40)53 + 

/ 80 40 44 \ 

(24-8N_ + 16N+2-32N+)5i,_2+( -YN+ + yN+2 + 28-yN_J5i,i + 

(48 - 16N_ - 64N+ + 32N+2)5l2 + (64N+ + 16N_ - 48 - 32N+2)52,i + 
(8N_ + 32N+ - I6N+2 - 24)5i,i4 1 + 

'^/{(B-yN_ + fN,2-^N,)5.+ (-8 + ^N_-fN,2 + fN,)52 + 

(-84N_-^N,2.fN,)5,}, 



iSiN) = 



C/ 1 (^(16^2 - 18)N+ + (64^2 - 2)N_ - 48^2 + 20 - 32C2N_2] Si + (siN^ - 34 + 
2N+ ) 52 + (8N_ - 12 + 4N+)53 + (8N+ - 8)5i,i + f 16N+ - 32N_2 + 64N 



(B.ll) 



' ■- -^/-j ' \ — --r -/~i,i ' 1 --r ■— i. ' - ■-■ — 

48 ] 5i,2 + (48 - 64N_ + 32N_2 - 16N+)52,i + (24 - 8N+ + 16N_2 - 32N_)5i,i j l + 

C.Q { (f N_. ^ 1Z5n,, ^ Ii?N_ - 24 - il^N,) 5, + (f N,. + f N. + 24N_. - 

4 - 56N_ J 52 + (24N+ + 16N_ - 8 - 32N_2)53 + (32N_ - 16N_2 + 8N+ - 24)5i _2 + 
(8 - 8N+)5i,i + (48N_2 - 24N+ + 72 - 96N_)5i,2 + (32N_ - 16N_2 + 8N+ - 24)52,i + 
(-24 + 32N_-16N_2 + 8N+)5i,i,i| , 



iSiN) 



CpUf 



(B.12) 



,^ 328^^ 256^^ 184^^ 

32 + --N+2 - --N_ + --N_2 ■ 



544, 



N+ 5i + 



21 



(32 28 32 52 \ 1 

_N+2-8-yN_-yN_2 + yN+jS2 + 2 + (8N_-8N+)S3| + 

^^^^ + 64C2')n+-32^2N+2+(^ + 64^2')n_-96C2-^ 



C/<^85_3 + 



32^2N-2J5i + ( -yN+2-40 + yN_2 + yN+-yN_)52+('-16N+2- 



32N_2 - 72 + 16N_ + 80N+ J 53 + (32N+ + 32N_ - I6N+2 - 16N_2 - 48)5i _2 + 

(48 - 32N+ + I6N+2 - 32N_ + 16N_2)^i,2 + 

88 32 1 

(48 - 32N+ + I6N+2 - 32N_ + 16N_2)52,i + y ^2 - y [ + 

C.«/{(fN-fN.-|-|N..fN,).,. 

The O (tts) coefficient functions in A'^-space read: 

41{N) =C^|(-3N+ + 6)Si + (4N_ + 4N+)52-9 + (2N_+2N+)5i,i| , (B.13) 



4'JW=Cf(2-2N_)Si, 

4!qW -4iiN) = C^(2N_ +2N+ -4)5i 



(B.14) 
(B.15) 



4^](A^) =q|(-16N- + 8N_2 + 8)5i + (16N_2-8N+ + 24-32N_)52 (B.16) 



+(8N_2 - 4N+ + 12 - 16N_)5i,i 



c[^l{N) = Q(16N_-8N_2-8)5i . 



(B.17) 



Finally, for the (a^) coefficient functions we give the even Mellin moments of c) J^^, c] g, Cj 



with / = r,L in Eqs. (lBJ8t - (lR23t and the odd Mellin moments of cV' in Eq. (EM 



41w 



,ps' /,g 



(B.18) 



5(A^-2) 



CpC^ 



e{N-3) 



c. 



560^ 704 



15 



^3 + 



145517 
810 



+ 



52^ 316^ 1273491 ^ 



810 J 



2354 „^ 16^ 

h8C3H C2 

81 ^ ^^^ 9 ^^ 



+ 



331 



C^^<J59^2-12C3+V25_4-16S_2 + 16S_3,i+48S_2,-2 + -^ + 



118 



+ 



22 



I2C2) N+ - y N+2 + (4^2 - ^"j N_ - 26 + y N_2) S2+( ls^3 - y - I2C2) N+ + 



4_8 



48 _ 279^ 



'48 48 279 \ /48 48 

^N+2 + (-63 - 12^2 + 8C3)N- - 48^3 + y N_2 + y^ ) ^1 + ( y N+3 - 36 - y N+2 - 



24N+ 

C/7C4 < 85-2 + 



30N+ j 53 - (22N_ + 138N+)54 - (16N_ + 16N+)Si _2,i + (32 - 48N_ - 48N+)5i,i,_2 H 

/ 48 48 48 48 

(_16 - 16N_ - 16N+)5i,_3 + ( - 8N_ - yN_2 + 16 - yN+2 + yN+3 + yN_3 - 

8N+yi,_2 + (56-28N+-4N_)5i,2+('64-(32C2 + 51)N+-(32^2 + 40)N_')5i,i + 

(-28N_ + 16-28N+)5i,3 + (-16-32N+ + 32N_)52-2 + (-16 + 24N_+4N+)52,i + 
(36N_ + 36N+)52,2 + (16N_ - 24N+)53,i + (8 + 20N_ + 8N+)5i,i,i + (24N_ + 

"" )5i,i,2 + (16N+ + 16N_)5i,2,i + (16N+ + 24N_)52,i,i + (24N_+24N+)5i,i,i,i| + 
((lf-3.C3-fc.)N.-^N..(li5_3.,_|,)._. 

.4C3-|N.-f)s,.((li^-S,)...|N..(f-S,)._-|..- 

38\ /23 24 24 H \ 

yJ52+(yN+ + yN+2-yN+3-4+yN_J53 + 36N+54-85_3,i-245_2,-2 + 

/ 24 24 24 24 

(8 + 8N_ + 8N+)5i,_3 + f 4N+ - yN_3 - yN+3 - 8 + yN+2 + yN-2 + 

4N_')5i,-2+((^ + 8C2V+-^+('^ + 8C2)n_')Su + ('i2N_ + 12N+- 

S\ 5i3 + (8 - 16N_ + 16N+)S2,-2 + (y N- + y N+'j ^2,1 + (4N_ + 4N+)S2,2 + 

/22 

(_4N_ + 4N+)53,i + (8N_ + 8N+)5i,-2,i + (-16 + 24N+ + 24N_)5i,i,-2 + f y N- 



— N_ + 
54651 



(_4N_ + 4N+)53,i + (8N_ + 8N+)5i _ 

yN+')5i,i,i-(4N-+4N+)Si,i,2 + (4N-+4N+)5i,2,i -365-4 + 7^2 + 10^3-^1 + 

/ 4 4 \ 457 " 

(-3N+--N_)5m,-2C2 + ^ + 8C3 



C, 



c?Un) = 



(B.19) 



d{N-2)Cpnf 



2462 16 



+ 



23 



e(A^-3)C^n 



118 4^ \ ^ /512 16. \ ^ /982 4^ \ ^ /16^ 



80\^ 188 „^ \ //^ 184\^ 128^^ 

2:yJN-2- — -8C2J5i+n8^2 + -^JN+ + — N+2 + 



184 



8^2 N_+32- 



y N^z) S2+(- y N_2 - y N+ + y N+2 + ^N^") ^3 + (-44N+ + 44N_)54 + 

^N,_32-^N.2-^N,2 + ^N_)5,_2+(f-28N_-^N_2-12N,- 

32 \ / 16 32 44 44 \ 

yN+2 J5i.i + (16 - yN+2 - yN_2 - yN+ + yN- J52,i + (24N_ -24N+)53,i + 

- y N+2 + 8 - y N_2 + ^N- + ^N+'j Si,i,i + (-8N+ + 8N-)52,i,i 



4'eW 



(B.20) 



5(A^-2) 

0(A^-3) 
316 



,,1408^ 1120^ 1406571 ri04^ ^„^ 26431 



15 



405 J 



81 



+ 



16, 



C/<; ( (36^2 + 34-16^3)N+-yN+2 + (184 + 32^3 + 208C2)N- 



^ + 96^2 + 64^3) N-2 - 148^2 + 48^3 ^ ^) ^1 + ( (^ + ^4^2) N+ + y N+2 + 

,..^ 274 \^ /144 
112^2 + ^ 1N_+' 



V 5 



64^2 ) N_2 - 72^2 - 136 j S2 + ( 166N+ - 230 + y N+2 



16, 



N+3 + 64N_ 53 + (- 132 + 88N_ + 44N+)S4 + (96 - 96N_ + 32N_2 - 32N+)5i _3 + 



64N_ -yN+3- 160 -yN_2 + yN_3 + 96N+ + yN+2JSi,-2+( (16^2 + 24)N+ + 
(-8 + 64^2)N_ + (8 - 32^2)N_2 - 24 - 48^2 j ^1,1 + (-16 + 16N+)5i,2 + ( - 96N_2 + 
64N+ - 192 + 224N_ j 5i,3 + (96N_2 - 12N+ + 92 - 176N_)52,i + ( - 128N_ + 

128 j 52-2 + (64N_2 - 128N_ + 96 - 32N+)52,2 + (-8N+ - 80N_ + 64N_2 + 24)53,i + 

(-64N_2 + 64N+ - 192 + 192N_)5i,i _2 + (-28N+ - 208N_ +96N_2 + 140)Si ,1,1 + 
(-32N_ + 16N_2 - 8N+ + 24)5i ,1,2 + (72 - 24N+ - 96N_ +48N_2)5i,2,i + 

(120 + 96N_2-40N+-176N_)52,i,i + (80N_2-40N+ + 120-160N_)5i,i,i,i| + 

8^\xT f 928 32^ \^ /440. ^^^„ 

-C2jN++(^- — -yC2jN+2+(^ — C2 + 24OC3 + 

^ 4438\_ _„ __ ._ \„ // 260 



Cy^Cf 



1934 



fV2C3 + ^ 



5410\^, / 248 



C2-96C3-^')n_2 + 70-216C3-56C2Vi +('('- 



24 



48^2 
/464 



N+ - ^N+2 + (92 - 64C2)N- + (^^+ 64^2 )n_2- 200 + 48^2)^2 + 

O^O /I /I oo \ 



■N_2- 



44 32 \ 

•N_ + — N+ - — N+2 - 36 J 53 + (248N+ - 320N_2 + 144N_ - 72)^4 + 



( 48N+ - 144 + 176N_ - 80N_2 j^l-s + [^^+2 - 



320^^ 80^^ 200^^ 

— N_ + yN_2--^N+ + 



48N,-144 + .,_ .._^_.,^3.,,. 3.,-, 3_ 3.,., 

,)5i,_2+((32C2-f)N+ + yN+2+(l28C2-^)N_+(^-64C2)N_2 + 

232 \ / 

— - 96^2 J 5i,i + (72 - 128N_ + 64N_2 - 8N+)5i,2 + f 96 - 32N+ - 1 12N_ + 

48N_2 ) Sl3 + (128N- - 64N_2 + 32N+ - 96)52,-2 + (^^+2 + ^^-2 + Y^+ " 

784 \ / 

— N_ + 80 J 52,1 + (96N_2 - 48N+ + 144 - 192N_)52,2 + ( 32 - 96N_2 - 32N_ + 

96N+ j 53,1 + ( 64N_ - 48 - 32N_2 + 16N+ j 5i,_2,i + (-32N_ - 16N+ +48)5i,i,_2 + 

/344 80 632 32 \ / 

f— N-2-yN+-— N- + 112 + yN+2j5l,l,l + (64N_2-128N-+96- 
32N+') 5i,i,2 + (32N_2 - 16N+ - 64N_ +48)5i,2,i + 



(24N+ + 24 - 32N_ - 16N_2)^2,i,i + (24 - 32N_ + 16N_2 - 8N+)5i, 1,1,1 



c?LiN) 



(B.21) 



5(A^-2) 



e(A^-3) 



C 



2 [48. 33 , ^ ^ 



24 2211 11 



10 



C/|r-10N+ + yN+2 + (-48C3 + 39)N_ 



+ 



f-fN_2 + 48C3)5i + 



^N+ - |n+2 + f N- - y N_2 - 2) 52 + (l2N_ + 32N+ - |n+2 + f N+3 - 
44 J 53 + (16 - 16N_)5i,_3 + ( 16N_ + y N_2 + 32N+ - y N+2 + y N+3 - y N-3 - 
48 J5i,_2 + (14N_ - 10 - 4N+)5i,i + (16 - 16N_)5i,2 + (16N_ - 16)5i,3 + 
(32 - 32N_)52,-2 + (4 - 4N_)52,i + (-32 + 32N_)5i,i,_2 + (-8N_ + 8)5i,i,i 1 + 



24 



lO^T 16^, / ,^ 389 \, 2023 ,,^ 24^^ \ /34 



24 16 218 \ / 16 16 

^N_2-yN+ + yN+2-^N_J52+(l6 + yN+2-yN+3-16N+)53 + 



25 



(8N--8)Si,-3 + 



T^-^- 5 



-N+3 + y N+2 + 24 - 8N_ + y N_3 - leN+'j 5i,_2 + 
(^ - y N-") 5i,i + (-8N_ + 8)5i3 + (16N_ - 16)52,-2 + (16 - 16N_)5i,i,_2} + 



CfHf 



4i(A^) 



fN-f4N0^-(l-l)^-G--l^... 



(B.22) 



Oft 

-—d{N-2)Cpnf + B{N-3)Cpnf 



f- 



^ + 8N_2 + ^N+-8N+2)5i + 



(^y N+ - ^N+2 - y N_ + y N_2) 52 + (24 - 24N_)53 + 
^N+ - y N_ - ^N+2 - 8 + Y^-i) Si,i + (8 - 8N_)52,i 



4'kA^) 



(B.23) 



5(A^-2) 

e(A^-3) 
104 



C. 



c 



96 1081 272 

yC3-— +— C^Q 



+ 



^"?N_ 



176 96 32 \ /64 

_ + _N_,--N,.-24N_ 5,+ -N- 



40 

T 



96 32 \ / 32 32 \ 

N_ + y N_2 + Y^N+2 J 52 + (48 - y N+3 - 48N- + y N+2 J ^3 + 

/32^, 32^^ 64^^ 64^^ 32^^ 32\ ,^^ ..^. ^. 

( yN+2 + y N- + yN_2 - yN-3 - yN+3 - y J 5i,-2 + (16 - 16N_)52,i + 

(8N+ - 56N_ + 32N_2 + 16)5i,i | + 

C,Q|f-64N+ + fN+2-^+256N_-^N_2)si + (fN+2-^N+-80 + 



688 



N_ 



352 \ 

— N_2 ] ^2 + (-96 + 128N_2 - 32N_)53 + (-64N_ + 32N_2 + 32)5i,_2 + 



(-64 + 128N_ - 64N_2)5i,2 + (64N_ - 64)52,i + ( y^N_ - 128 + yN+2 

464 64 \ 

— N_2 - y N+ 1 5i,i + (-32N_2 + 64N_ - 32)5i,ia 



J2) 



CYLJN) 



(B.24) 



5(A^-1)C^ 
0(iV-2) 



11 



— Q - -rif 



2 
3^ 



12)^3 + 
C/| ((y -40^3- 12^2^+ + (-16-40^3 -12C2)N- 



19 + 48^3 5i + 



26 



((12^2 -92)N+ +78 + (4^2 -92)N_)52 - (94N_ + 66N+)S4+ f 20N_ -90N+ + 
16N+2 - 12 j 53 + (24 + 12N_ - 12N+)5i,2 + (-12N_ - 12N+ - 16)5i3 + ( 16 - 
32N_ j S2 -2 + (-8N+ + 44N_ - 24)52,i + (36N_ + 36N+)52,2 + { - 32N_ + 16 - 
32N+ j 5i,_3 + (16N_2 - 24N+ + I6N+2 - 24N_ + 16)5i,_2 + ( (-32^2 - 33)N+ + 

28 + (-22 - 32C2)N_ j 5i,i + (24N_ + 16N+)52,i,i - (16N_ + 16N+)5i _2j + 

(-32 - 16N_ - 16N+)5i,i -2 + (28N_ + 16N+ - 8)5i ,1,1 + (24N+ + 24N_)5i,i,2 + 
(16N+ + 16N_)5i,2,i + (24N+ + 24N_)5i, 1,1,1 + 725_4 - 165_2 + 165_3,i + 

485_2,-2 - 8N+53,i + ^ - 12^3 + 59^2! + 



c^c 



F*-A 



193 



(-yC2 + ^-12C3JN,+ ^--C2 + ^-12C3JN_---24C3j5i + 
■8C2)n+-^+(^-8C2)n_)52+(-4-^N_ + ^N+-8N+2)s3 + 

/ 9 Afs \ 



(8 + 4N_ +4N+)5i,3 + (16N_ - 8)52,-2 + ( " ^N_ + y N+ j ^2,1 + (4N+ + 4N_)52,2 + 
36N_54 + (4N_ - 4N+)53,i - 85„3,i - 245„2 -2 + (16N_ + 16N+ - 8)5i,_3 + ( 12N_ + 

12N+-8N_2-8N+2-8)si,_2+('('8C2 + ^')n+-36+('8C2 + ^)n_')5i,i + 

/22 22 \ 

(8N_ + 8N+)Si,_2,i + (16 + 8N_ + 8N+)Si,i,_2 + ( y N+ + ^ N_ j Si,i,i - 365_4 + 

+ (-4N+-4N-)5i,i,2 + (4N++4N-)5i,2,i-^ + 10^3+7^2! + 



85-2 



CpUf 



'^r 131\ ^ 62 (A^ 

3C2-^)N+ + -+(-C2 



302\ \ /16 22 26 \ 



2 2 

- -N+ - -N 
3^3 



^3+(-^N+-^N_)52,i+(-yN+ + 4-^N_)5i,i + 
(-lN,-lN_)5i,i,i + f + 8C3-2C2[. 

Appendix C: The exact x-space results 

Here we write down the full x-space results of the time-like splitting functions and the coefficient 
functions cj-, c^ and ca up to second order in a^ = as/(47r), expressed in terms of harmonic poly- 



27 



logarithms in the notation //„7i,...,m,^,(x), nij = 0, ±1 of Ref. [67]. Below we use the short-hand 
notation 

'^0,...,0,±1,0,...,0,±1,...W = ■^±(m+l),±(;7+l),...W: (C-1) 

m n 

suppress the argument x for brevity and define 

/JqqW = 2{\-x)-^-\-X, (C.2) 

Pgg{x) = (l—x)^ +x^ —2+x — x . 

All divergences for x ^ 1 are understood in the sense of +-distributions. 
The well known results at LO are: 

<^ W = Cp (2pqqW +36(1 -x)) , (C.3) 

PJg^ (jc) =nf{2-Ax + Ax^) , (C.4) 

Pf^ [x) =Cp(-A + ^ + 2x\ , (C.5) 

^i? W = Q (4;'ggW + y5(l -x)^ -nf\h{\ -x) . (C.6) 

The NLO splitting functions read"^: 

Pis^+W = (C.7) 

C^'ji^qq W(-8H2 + 6Ho - 16Ho,o - 8Hi,o) +7?qq(-x) (-8C2 + 8Ho,o - 16H_i,o) + 

(-6 + 2x)Ho + (4 + 4x)Ho,o-4 + 4x+('-12C2 + ^+24C3')5(l-x)| + 

r / 134 22 \ 

QQJPqqW (4Ho,o-4C2 + -^ + yHoj +Pqq(-x)(-4Ho,0 + 8H_i,o+4^2) + 

|(l-x)+(yC2 + y-12C3)6(l-x)} + 
«/C^{i^qqW(-^Ho-^)-^(l-x)+(-^-^C2)5(l-x)|, 

Pi^'{x) = (C.8) 

C/|pqqW(6Ho-8H2-8Hi,o-16Ho,o)+Pqq(-^)(16H_i,o-8Ho,o + 8C2) + 

(-14x-22)Ho + (4 + 4jc)Ho,o + 36jc-36+Q- 12^2 + 24^3 j 5(1 -jc)| + 



'^We quote f qg and f gg here for completeness (see discussion in Section|3i. 
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f / 22 134\ 

QQJ/'qqW ( -4C2 + 4Ho,o + yHo + — J +(-4C2-8H_i,o + 4Ho,o)Pqq(-x) + 

^ f / ^/^ 4^^ 20\ 8 8/1 8^ \ ^,^ , 
CpHfip^dix) I -:^Ho-— ) +^x--- [ ^ + :^(:,2 )o(l-xJ 



p!^\^\x) =nfCp{S{\+x)no,o 



9/3 3 

32 



3 3 



vC^ + 20 + 36x |Ho + 



16(x-1)(14jc2 + 23jc + 5) 
9x 



(C.9) 



Pi\^) 



HfCp I (8x - 10 + 8x2)Ho + (8^2 + 12 - 8x)Hi + (8x - 4 - 16x2)Ho,o + 



46;c+16;c^ + 8)(H2-3Hi,o-Hij-^2)-40;c^-24 + 46;c^ + 



^ , , 136 16 2 16, „ ,40 40 . ^^^.. 
nfC, <! ( — —X - —x^ - y )Ho + {—X - —x^ - y )Hi + 



nf^A 



(-16x+16jc2 + 8)(Hij-2H2 + 2Hi,o) + (-16jc-8-16jc2)H_i,o + 

.o .o N ..^ 4(569jc2- 20+ 178^3 + 13x)l 
(8 + 48x)Ho,o-16C2X+^ — '-i + 



16 2 8 16 

y^+3-y 



2 J / ^V + ^-^x)(Hi-Ho) 



9x 

40 32 32 2 

1 X X 

9 9 9 



(CIO) 



piW^) 



C^2J (_32 + 2jc)Ho-8jcHi + ^^^:^^^^^^-t^(Hi,i +2Hi,o-2H2) + 



8-4x)Ho,o+18;c-2^ + 

4(27;c2^8;c3^ 24^-18) 



Cj^Cf 



(C.ll) 



3x 



Ho + 8xHi + ^^ ^^+^' + ^) (H2-3Hi,o-Hi,i) + 



8(2;c + x2 + 2)^^ S{2x + 3x^ + 4)^^ ^^. 4(-45x + 9x2 - 17 + 44jc3) 1 
^ ^-i'« ^ Ho,o + 16C2 ^^ j 



P^'J(x) 



rifCp 



4(21;c2 + 8;c3 + 15;c + 8) 



3x 



Ho + (8x + 8)Ho,o- 



8(x-1)(41j«;2 + 14x + 23) 



9x 



Q i/'ggW 



268 



25(1 -x)^ + 



88, 



^ I6H2 - 8^2 + yHo - 24Ho,o - 16Hi,o ) + 



(C.12) 



29 



4f44— ]]x-\-2'^x^) 
(-8^2 + 8Ho,o - 16H_i,o)pgg(-x) - ^ 3^^ ^ Hq + (-32 - 32x)Ho,o + 

2(.-l)(134.^-109.+ 134) /^^^^^^X'^^^_^1^ 



9x V 3 ^ 

^r/ 16^^ 40\ ,, /8 8 \ ^ 4(x-1)(13jc2 + 4jc+13) 8^,^ , 

The coefficient functions at order C'(as) read: 

c'^^'l(x)=cJp^^ix)(-^+mo-2Ui]-^x+^ + {Sl^2-9M\-x)] , (C.13) 

4;jW = 2C^ , (C.14) 

4'IW-4!!,W=Cf(2^-2), (C.15) 

c?J.) ^ C, f8(-2. + 2+.^)^^ _ 4(-2. + 2+.^)^^ ^ 8(-l+^A ^ ^^_^^^ 



c« W = C,^ . (C.17) 

The coefficient functions at order O (a^ ) read: 
41W = (C.18) 

C/|pqqWU20C2-53)Ho + 6^2 + 6H2 + yHi+4H3 + 12Hi,o + ^-18Hi,i + 
33Ho,o - 76^3 + 20Hi,o,o + 24Hi,i,o - 24Hi j,i + 16Hi,2 - 36H2,o + 24H_2,o + 20H2,i - 
80Ho,o,o j +Pqq{-x) ( - 32H_i^2 + (-8C2 - 8)Ho - 8H_2,o - 8H3 + 36Ho,o,o - 

32H_i,_i,o-24H_i^o,o + 16H_i^2 + 28^3V(-8-^-8)C2H-i + f(8x + 8)C2 + 



-59x2 + 48 + 48x3 + 71jc\„ /, „ „ ,. 75 53\^^ ,^„ ^ ,„ 

)Ho+ (-8 + 8x)C2-^^+^ Hi + (18 + 2x)H2- 



5jc / V '"2 2 

(40x-16)Hi,o + 16H_2,o + (-10jc-2)Hi,i + (4x + 4)H2j + (-16-16x)H_l_lo + 
(8x + 8)H_Lo,o + (22x + 22)Ho,o,o + (-8x + 8)Hi,o,o+(6x-yx3_igV2 + 

3(271jc-101jc2-64 + 64x3) , /331 , ^^,2 



(8 + 24x)C3 + ^ 20^ ^ ^+^— + 30^2^-78^3 -39C2J 5(1 -x)j + 

C^Q{/'qq(^)((^-16C2)Ho-(4C2 + ^)Hi-4H2,o-^-fHu + 14C3- 
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11 22 \ / 

12H_2,o-yHo,o-8Hi,o,o-4Hia,o + yH2 + 18Ho,o,o + 4Hi,2]+Pqq(-x)(l6H_iC2 + 

(4^2 +4)Ho - 18Ho,o,o + 4H_2,o - 14^3 + 12H_i,o,o - 8H_l2 + 16H_k_i,o + 4H3 j + 

,, .^. TT 2(-13jc2 + 36-108x + 36x3)^^ {.^ ^ .^ 139 85 \^ 
(4jc + 4)^2H_i - ^ — ^Ho + ( (4-4^)^2 - -g- + y-^ ) Hi - 

on-2,0H p-2 n-1,0- T^l-J+OX Jhlo,0 + 

(8x + 8)H_i,_i,o + (-4-4x)H_i,o,o + (4x-4)Hi,o,o + ^^(-5 + 6x2)^2 -8^3^- 
(2039x2 -864 + 864x3+ 441 Ix) / 49 5455 ^40 215^ \^, ,1 



r / 16 4 29 

CpnAp^^{x) ( - yHo + -Hij + —Hi + 



247 4 2 \ / 2 



5 M„ 43 71 M, 457 38, , .,, . 
3 + 3. H,+- + -.+ -C, + — --t|8(l-.) 



4'L(^) 



C,n,[ ((,6+ ,6.)fc - 8('^^ + 7^;;^ + 39.^) ^ ^^ _ ^^^^^^^^^ _ 

4(x-l)(8x2 + 25x + 2) 4(x+l)(4x2 + llx-8) 16(x+l)3 



3x 

2(32 -21x2 + 8x3 -21x) 
3x 



H, 



0,0 



3x 

4(x-l)(4x2 + 7x + 4) 



3x 



3x 
Hij + (8 + 8x)H2,i + 



(44x + 44)Ho,o,o + (-8-32x)C2 + (16 + 16x)C3 + 



2(x- 1) (256x2 + 571X + 40) 
27x 



(C.19) 



4^iw 



C/|^^^^^C2H-i + (^(-16x + 32)C2 



2(-209x- 72 +131x2 + 8x3) 



5x 



Ho + 



(C.20) 



8(5x2-10x + 6)^ 8(3x2-l)\ 4(3x2-20x + 24) 8(-2x + x2 + 8) 
■^ -(,2 + ^ ) Hi + ^ ^H2 - ^ ^Hg - 



128H_2,o + 



16(-20x2-30x3+x5-4) 



5x2 



H-10+ ( -64 + 166X-— x^ )Hoo-16Hi^o^ 



4(-28x + 7x2 + 24)^ 24(-2x + 2 + x2)^ 32(-2x + 2+x2)^ 
Hij H Hi^2 H2,o + 

8(-10x + 12 + 5x2) , 64(x+l)2 32(x+l£ ...^««m ^ 

-H2,i H H_i -1,0 H_i^o,0 + (-44x + 88jHo,o,o + 



X 



32(2x2 -4x + 3) 



H1.0.0 + 



8(-2x + 2+x2^ 



H, 



40(-2x + 2+x2) 



X 



X 



H1.1.1 + 



31 



16 3 n^\y 16(4 + 5x2^2x) 2(302x + 8x3 -77^2 -158) 

5 J X 5x 

^^\[ 16(8x + Jc2-2). 4(579x + 372 + 64x3 + 129jc2)\ 

c^Q|(^ — -, — ^2+ ¥x n^ 

8(4 + jc2 + 2jc)„ ^^ /8(3jc2-6x + 4)^ 4(8;^^ +^2 - 89 + 59x)\ ^ 

-Q2H-1 + C,2 H ^ Hi - 



+ 



X \ X 3x 

16(3x2 + 18jc - 31 +2jc3)^^ 32(4jc + 3 + 3^2)^^ 32(2x + 2 + x2)^^ 
z H2 H H3 H H_2,o + 

8(4jc3 + 10 + 30x + 21jc2) 16(2jc + 2 + x2) 4(24jc-3jc2 + 116 + 8x3) 
^ tl_i H_i 2 Ho,o + 

8(8+x2-8x)^ 8(36x-6x2+4x3-43)^ 16(-2x + 2 + x2) 
Hi^o H ^ Hij H Hi^2 — 

lA. JT^ ■A' 

48(-2x + 2+x2) 8(3x2 + 6x + 2) 

H2,o H2,i + (-16x-32jH_i -1,0 + 

16(5 + 6x + 3x2) 8(40 + 22x + 31x2) 16(2x2-4x + 3)^ 
ti_i,o,o Ho,o,o tii^o,o + 

32(-2x + 2 + x2)^^ 8(-2x + 2+x2) 8(-7 + 8x2 + 4x) „ 
Hi, 1^0 Hi.i J H 1,2 — 

8(17x2 + 30 -8x) 2(486x- 2219 + 464x3 + 1431x2) 

X 27x 

41w = (C.21) 

'"^^-^ + ^^^7^ + '"-''' H_i.o+('^x3- 12 + 32x)ho,o-16H|.o + 8Hu-32H_,,o + 
32H^,,_,,„ - 16H_,,„,„ - ,6H,,„,„ + (4 - f .' - 32.) C. + 48- 147.- 18x^ + 32.^ | ^ 

C,c,{ - 8H_,C. - ?tl^^^^i^t2!f±^H„ + (f - 8fc) H, + 

8(2x^+5x2-3 + 10x3) 16 , o , 

^ —^ ^H_i,o - yx(x2 + 5)Ho,o + 16H_2,o - 16H_i,_i,o + 8H_i,o,o + 

1h,-^ + 1h„+1. 

3 9 3 "3 



Cpnf{ - -Hi - -^ + -Ho + -X [> , 



'L,psW ~ 

[ 8(-6x2+x3-6x + 4) ^ 8(x-l)(x2-2x-2) ^ 
C^n^ 3^^ Ho Hi - 8H2 + 



32 



(C.22) 



24Ho,o + 8C2- 



8(jc-1)(3jc2-10jc-3) 



4'^(-) 



C, 



(C.23) 



8(;c + 4)(;c-l)^^ 8(28x2 -36 + 3jc + 4x3) 



-H 



15x 



32(x+l)(x^-x3+x^-6x + 6) 



C^Q 



32(x+l) 



15x2 
16(22- 6x2 -21x + x3 



H-i,o+( 16 + ^x3)^2 



Ho-(^x3-48JHo,o-16H2 + 

8(x-l)(4x2-27x-36) 



H 



-1,0 



3x 

32(3x + 4) 



-Ho + 



32 
15 
16(x-l)(x2-2x-29) 



3x 



15x 
Hi+64H2 + 



32(2x-l)„ 32(x-l)(x2-4x-14) 

Q2H ^ 

X 3x 



64(x-l) 32(x-l) 

nO,0 H ni,0 ni J 



X 



+ 



^Sw-4'iw 



(C.24) 



C/|pqq(-^)('(16C2 + 16)Ho + 48H_i,o,o-32H_i,2+16H3 



72Ho,o,o + 16H_2,o + 64H_i,_i,o - 56^3 + 64H_i^2 



2(153x + 24 -107x2 + 24x3) 



5x 



Ho + 



((16 - 16x)^2 + 18x- 18)Hi + (20+ 12x)H2 + (16 - 16x)Hi,o + (-8 + 8x)Hi,i + 

( l^-H^.m , 16(x+l)(3x2+x + 3)(x2-3x+l) 

(-16 + 16xjH_2,o H ri H_i,o + 

48 \ / 48 

-20-44x+16x2-— x3 Ho,o-(16-16x)Hi,o,o+ -16x2 -20x-28 + — x^ ) ^2 + 



;,6-,6.)C3-(i^il(l?4^H^i±l^,+ 



CfQ|p,<,(-i)((-8C2-8)Ho + 28i;3 + 36Ho,o,o-8H3-32H_i,_,,o + 16H_i,2- 
24H_,„,„ - 8H_,„ - 32H_,C.l + 4(18-104.^ + 18.^+161.,) ^^ /^_^ ^ ^^^^^ , 



15x 
|.^^)h.+(-8.-8)H, + (8-8.)H_,„- ^'-+')P-'^+;;/)'-'-^-^" h_,o + 

16x+ 16 - 8x2 ^ ^^3 j fj^^ ^ ^g _ 8;,)Hi + /^ 8x2 ^ ^^x^ 12 - yx^ j ^2 + 
(-8 + 8.)fa+'^-'"^'^^:.^'"^ + ^'^'K 



C^Hy 



45x 

4 4 \ ^ /4 4 \ ^ 38 38 , 

h-x Ho+ X Hi xH \ . 

3 3 / " V3 3 / 9 9^ 
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